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^O ■ Abstract 

(D ; 

pL^ ! We show that massive sterile neutrinos mixed with the ordinary ones may be 

j — ■ produced in the early universe in the right amount to be natural warm dark matter 

particles. Their mass should be below 40 keV and the corresponding mixing angles 

^^ | sin 2 2# > 10~ n for mixing with v^ or v T , while mixing with u e is slightly stronger 

£> ' bounded with mass less than 30 keV. 

en ! 

oo ■ 

§ ■ 1 Introduction 

o: 
o- 

There seems to be convincing experimental evidence for non-zero neutrino masses and 

mixing angles (for a review see e.g. |l|), and if all the present day data are correct, there 

Ph. must exist at least one sterile neutrino species. These neutrinos should be very light 

(sub eV range) and hence contribute negligibly to the cosmological energy density 
> ■ 

'$ : n„/i a = ^- , (i) 

^ . 92eV V ; 

if they were produced with the equilibrium number density in the early universe at a 
temperature below ~ 10 MeV. If their number density was smaller than the equilibrium 
one, then the permitted value of the mass could be correspondingly higher. 

One could easily envisage more than one sterile neutrino species. The masses and 
mixing angles of these extra neutrinos are essentially free parameters. If we consider 
sterile neutrinos with masses 10 — 200 MeV, then big bang nucleosynthesis and energy 
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loss arguments for SN 1987A allow one to exclude mixing angles in the range sin 2 2# = 
lCT 1 — 10~ 12 H, where both the upper and lower limits vary as functions of mass. On 
the other hand, direct terrestrial experiments exclude supplementary mixing angles in 
the range sin 2 2# = 0.001 — 1 ||], and the lower limit is weaker for small masses, m ~ 10 
MeV. We should mention that the excluded regions only overlap for some masses, and 
there is still some non-excluded parameter space for m < 40 MeV and sin 2 2# = 0.01 — 1. 
The hypothesis that sterile neutrinos could make a considerable contribution to cos- 
mological dark matter has a rather long history. The idea that right-handed sterile 
neutrinos may form warm dark matter was briefly discussed in ref. || and was further 
pursued in the paper || . In more detail warm dark matter cosmology was considered in 
ref. . There are some other warm dark matter candidates discussed in the papers || || . 
Sterile neutrinos coming from a mirror world may also make WDM in our universe; they 



were discussed in the papers ||10|| . More models and references can be found in the re- 



cent works [11, O]. A dark matter model with sterile neutrinos but with a non-thermal 



spectrum was considered in ref. [|13|] . Such neutrinos could be produced by the resonance 



oscillations in the early universe in the presence of a large lepton asymmetry. This model 



was further considered in ref. |14j], where constraints originating from consideration of 
decays of v s , especially of the radiative one, were presented. A detailed analysis of the 
essential physics of this model, namely the importance of the mixing angle suppression 



in the early universe, is made in ref. |L5] 



In this paper we find the allowed values of mass and mixing angle of a sterile neutrino, 
z/g, so that the latter could be a dominant dark matter particle. We will consider the 
production of v s in the early universe, calculate their energy spectrum, discuss their 
different decay modes, and derive bounds on mass and mixing angle from cosmology 
and astrophysics. 



2 Production of v s in the early universe 

Let us consider for simplicity a two-neutrino mixing scheme, where one of the active 
neutrinos, v a = v e , v^ or u T , mixes with a heavy mainly sterile neutrino, u s , 

v a = cos 9 v\ + sin 9 v 2 , 

v s = — sin 9 V\ + cos 9 v 2 > (2) 

where v\ and v-i are assumed to be the light and heavy mass eigenstates respectively, 
and 9 is the vacuum mixing angle. We will consider small mixing angles, and hence 
sometimes refer to the light neutrino mass eigenstate as the active neutrino and the 
heavy one as sterile neutrino. 

We assume that sterile neutrinos were initially absent in the primeval plasma and 
were produced through the mixing with active ones. As will be clear shortly, the 
temperature suppression of the effective mixing angle changes the production rate at 
high temperatures from T 3 to T™ 9 , so with small mixing angles the sterile neutrinos 



would never have reached equilibrium. As was estimated in ref. [[L6[] the maximum ratio 
of the production rate to the expansion rate takes place at rather low temperatures, 
T « 0.1GeV(m/MeV) 1//3 . The production rate is usually approximated as |16|, I7| 



where H is the Hubble expansion parameter, T is the plasma temperature, and T\y is 
the decoupling temperature of the active neutrinos, which is approximately taken about 
TV — 3 MeV. Instead of this approximate equation, below we will write down and solve 
the exact momentum dependent Boltzmann equation, taking into account the processes 
of production of v s but neglecting inverse reactions. The latter are not important if the 
number density of v s is small. In this more precise approach the question of the value of 
T\y never appears, it is solved automatically. Another advantage of our approach here 
is that it permits to calculate the energy spectrum of u s , while the previous method 



permitted only to estimate the total number density. Before doing these calculations it 
may be instructive to make the standard simplified estimates and later compare them 
with the exact results found below. 

The mixing angle sin 2 26* M is suppressed at large temperatures due to matter effects, 
and for v^ or v T mixing it can be written as (TT 



sin M ~ 1 + Q g x 1Q _ 19 ( T7 / MeV )6( 5m 2/ MeV 2)-i > y ) 

where the coefficient in front of the second term in the denominator was obtained by 
a rather arbitrary procedure of thermal averaging of the factor (E 2 ) w 12T 2 , entering 
the ratio of the neutrino refraction index to the vacuum term 5m 2 /2E 2 . We see from 
this expression that matter effects become essential and suppress the mixing for T 7 > 
0.15 GeV (m/keV) 1 / 3 (a similar argument was made in ref. || for the left-right neutrino 
mixing). For the [y e — zy s )-mixing the factor in the denominator should be 2.7 instead 
of 0.8. We will assume here that 5m 2 = m 2 , — m 2 is positive (specifically we assume 
mi <C 777.2), and if instead the active neutrino is heavy the analysis somewhat changes (see 
refs. 0,0). 

For the energy dependent calculations we need the expression for the matter effects 
in the denominator of eq. (|j) prior to averaging over the thermal bath. The latter can 
be read off from the relevant equations of refs. |TB, 113 



. sin2fl 

Sm M ~ 1 + 3.73 • 10- 20 c 2 m(MeV)" 2 (y 2 /x 6 ) ' ^ 

where the v s mass, m, is measured in MeV and we used the expansion parameter of 
the universe, a, to introduce the new variables, x = lMeV x a, y = Ea, and neglected 
a possible entropy release so that the temperature drops according to T = 1/a. The 
numerical coefficient c 2 depends upon the neutrino flavour: c 2 = 0.61 for v e and c 2 = 0.17 
for v T and v^. However, for the temperatures close to or above the muon mass c 2 becomes 
the same for v e and v^. 



The Boltzmann equation describing the evolution of the sterile neutrino distribution 
function, f s , in terms of these new variables takes the form 

xHd x f s = J co u , (6) 



where the collision integral is given by 

1 f d 3 p 2 



'coll 



\A\ 2 f 3 hdr 3A , (7) 



IE J (2vr) 3 2£ 2 

where dr 3A is the phase space element (together with the energy-momentum ^-function) 
of the particles l 3 and l± in whose collision v s is produced, 

h + h -> v s + h , (8) 

and f^ are their distribution functions. We assume that the latter are equal to their 
equilibrium values and then the conservation of energy gives f^f^ = exp(— yi — y 2 ) in 
the Boltzmann approximation. Integrating the probabilities of all the relevant processes 
over phase space (see the appendix) allows to find the collision integral and to solve 
equation @ analytically, giving 

/. = 3.6 • 10 8 sin 2 9 (l + g 2 L + g\) c~ 2 l/2 m(MeV) {^\ f a , (9) 

where f a is the distribution function of any of the active neutrinos, and g* is the number 
of relativistic degrees of freedom at the time when the sterile neutrinos were produced. 
Subsequent to the production there will be a dilution of the active neutrinos relative to 
the sterile ones. This is described by another factor (g*/10.75). 

The coefficient relating f s to f a in eq. (H) is independent on the energy of neutrinos, 
so the spectrum of v s remains the same as that of active neutrinos. This is somewhat 
surprising because the reaction rate is proportional to the neutrino energy. However, 
for smaller E the rate becomes efficient at higher temperature, as one can see from 
the expression (||) describing the suppression of neutrino mixing in matter. This effect 



compensates the factor of momentum, y, in the kinetic equation. However, at very large 
temperatures, T ~ Mw,z ~ 100 GeV, the weak reaction rate drops down, and hence the 
spectrum will be somewhat distorted at very small j/'s. 

It is interesting to compare the accurate results presented above with the simplified 
calculations based on the solution of the following approximate kinetic equation 

Hxd x f s = ^sm 2 29 M T w f a , (10) 

where the mixing angle and interaction rate can be taken from eqs. (Q, [|). This equation 
is easily integrated, and we find that the result for the total number density of v s agrees 
within a factor of 2 with the more accurate result (Q). 

Up to now we have seen how the produced amount of sterile neutrinos depends on 
the mass and mixing angle, so let us instead ask: how many sterile neutrinos should be 
produced in order for them to be a dark matter candidate? Let us take Qdm — 0.3, which 
means that we must demand p s = 3/i 2 keV/cm 3 . Using h = 0.65 and ?7,^ oday = 100/cm 3 
one finds 

„ /MeV\ /n DM \ ( h V 
«. = 127X10-V(— )(-g){^) • (ID 

and comparing eqs. (^) and flTT|) one obtains 

in i4 cj /2 /10.75\ 1/2 /MeV\ 2 . , 

sin 2 ^ = 3.6x 10" 14 - 1 ^ . 12 

(1 + 9l + 9r) V 9* J V m J 

This equation thus describes a line in mass-mixing parameter space, where the sterile 

neutrino must lie, if it indeed is the dominant dark matter particle. Let us now see how 

decay processes and supernovae can further restrict this parameter space. 

3 Decay 

The mixing couples the heavier v 2 to the Z-boson, and allows for the decay channel 

v 2 ^vi+Z + i, (13) 



where v\ is mostly an active flavour and I is any lepton with a mass smaller than half 
the mass of the heavy neutrino. This mixing angle can be translated into decay time 

r= 1Q5/(m) sec (14) 

m(MeV) 5 sin 2 2# ' { } 

where f(m) takes into account the open decay channels (for m < 1 MeV only the 

neutrino channels are open, and f(m) = 0.86, while for m s > 2m e the e + e~-channel is 

also open and f(m) = 1). Now, for the sterile neutrino to be a dark matter candidate 

we must demand that it does not decay on cosmic time scales, which means r > 4 x 10 17 

sec, and hence from eq. (|i"4|) we get 

sin 2 2£ < 2.5 x 10- 13 -^- . (15) 

m(JVleV) 

We can, however, get an even stronger bound by considering the radiative decay 

v s -> v a + 7 , (16) 

where v a is any of the active neutrinos. This decay will contribute with a distinct line 
into the diffuse photon background near m/2. The branching ratio for the reaction flTBD 
was found |B| to be: BR « 1/128. The flux of electromagnetic radiation from the 



decay was calculated in the papers p(| |2l[] (see also refs. [14], P2[). In the case of a large 
life-time, larger than the universe age, and of the matter dominated flat universe the 
intensity of the radiation in the frequency interval du is equal to 

dl = (BR) ^- ? ™ (17) 

1 ; Ht s (m s /2) 3 / 2 l ; 

where n^' is the present day number density of v s and H is the Hubble constant. Here 

we neglected corrections related to a possible dominance of the lambda-term in the latest 

history of the universe. 

In the energy range interesting to us a rather conservative upper limit for the flux 



can be read of the figure of ref. |23 



dT ni ZlMeV\ _ 2 _ x _! , 10 . 

— - < 0.1 — — - cm 2 sr ^ec L (18) 

d\ I \ E J 



and taking the values Q s = 0.3 and h = 0.65 we find: r > 4 x 10 22 , which together with 
eq. (|i~4"D leads to the bound 

sin 2 2# < 2.5 x lO" 18 f } m \^ . (19) 

m(MeV) 5 v ; 

A mass dependent bound can be found by considering the energy loss argument for 

SN 1987A, for m s < 3Tsn ~ 100 MeV. Sterile neutrinos produced due to mixing with the 

active ones inside the supernova would carry away too much energy, hence shortening the 

explosion. The excluded mixing angles have been calculated several times for SN 1987A, 

and the results are about sin 2 2# < 3 x 10~ 8 for v„ or v T mixing H] and sin 2 2# < 10~ 10 



for the v e mixing |24| for masses of the order MeV. For smaller masses, m < 40 keV, 

this bound weakens substantially due to matter effects, since the matter mixing angle 

can be expressed as [[Z5|] 

• 2~n sin 2 20 , . 

sm 2 26 M = . (20) 

sin 2 2fl + (cos2fl + 1.4 ■ 10 3 (^-)) 
On fig. 1 one sees, that this SN-bound is weaker than the diffuse gamma background 
limit for all neutrino masses. 

Plotting the equation describing the production, eq. (|12"|), together with the bounds 
from SN 1987A and radiative decay, makes it clear that (u a — z/ s )-mixing (see fig. 1) 
as the producer of dark matter demands that the neutrino mass must be smaller than 
40 keV for v^ or v r mixing (solid lines), and slightly below this values for v e mixing 
(dashed line), and with corresponding mixing angle sin 2 # > 10~ n . The two thinner 
full lines in the figure allow for a factor 2 uncertainty in the amount of dark matter, 
Q DM = 0.15 - 0.60. 

A lower bound on the dark matter particle can be obtained from large scale structure. 
A small mass of the dark matter particle will erase structure on small scale, and present 
day data seem to exclude masses smaller than a few tens of eV. 

Several comments are in order here. First we must check that the sterile neutrinos 
are indeed relativistic when produced. This is the case, because the temperature T max is 

8 
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Figure 1: Bounds from (u a — i/ s ) -mixing. The middle full line describes the mass-mixing 
relationship if sterile neutrinos are the dark matter for [v T — z/ s )-mixing. The two other 
full lines allow a factor 2 uncertainty in the amount of dark matter, Qdm = 0.15 — 0.6. 
The dashed line is for [y e — zA,)-mixing. The hatched region for big masses is excluded 
by the Diffuse Gamma Background. The region above the dotted line is excluded by the 
duration of SN 1987A for [y T — z/ s )-mixing. 



about 1.3 GeV for m = lMeV, and about 0.13 GeV for m = IkeV. Further, the dilution 
factor is somewhere between 1 and 4 depending upon whether the production happens 
before or after the QCD transition, and can thus enlarge the allowed region slightly 
compared to the figures, where we for simplicity used g* = 10.75. Looking at eq. (Q) it 
seems that the sterile neutrinos follow an equilibrium distribution function. This is not 
quite the case, because the small momentum neutrinos are produced first, and hence 
their relative importance is increased by the subsequent entropy release (which dilutes 
the active neutrinos). A different non-thermal effect can appear for (z/ M — z/ s )-mixing, 
since the factor C2 is 0.17 when the /i's are absent (for T <C m^), whereas it grows to 
C2 = 0.61 when the muons are fully present in the plasma. This means that bigger 
momenta will be produced with a factor 1.9 more efficiently 



9 



4 Discussion and conclusion 

The model considered in this paper, is undoubtly the simplest, oldest and, as we have 
seen, very natural for warm dark matter. The value of the v s mass can in the future 
be found from the detailed analysis of large scale structure formation. As we have seen 
this model permits masses about keV which may be interesting for the galaxy forma- 



tion problem |11|| , and furthermore recent N-body simulations of large scale structure 
compared with observations of the number of satellite galaxies have in fact indicated 
that the dark matter particle may have a mass about keV ||27|| . It is interesting to note, 



that a keV sterile neutrino has been suggested as a possible explanation of the observed 
pulsar velocities fSBR . Finally, better observations of the diffuse 7 background around 
keV energies should be able to cut away more of the parameter space, or potentially 
make an indirect observation of dark matter. 
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A Solving the Boltzmann equation 

All the relevant processes were presented in table 2 of ref. ||. There are two kinds of 
matrix elements, namely (P1P2XP3P4) and (piP4,)(p3P2), an d one finds from the integral 
over phase space that 

dTuip^iPaPi) = 3 / dr 34 (p 1 p 4 ){p3P2) = g (21) 



10 



Now one can count all the relevant processes, integrate over momenta and find 

5x2 4 



Hxd x f s 



3tt 3 



sin 2 9(l + gl + g 2 R ) G 2 F E x T A f a . 



(22) 



where Ex = 4.5 x KT 22 ($%) x' 1 MeV and G F = 1.1664 x 10~ 5 GeV. With the variable 
£ = y/x 3 and j3 defined in eq. ([5]) the suppression of mixing angle is 

sin 29 



sin29 



M 



i + p 2 e ' 



(23) 



and with the integral 



dti 



i y _4 i 



(24) 



we find the result in eq. 



References 



[1] J. Ellis, Nucl. Phys. Proc. Suppl. 91 (2000) 503 f|iep-ph/0008334 . 



[2] S.S. Gerstein, Ya.B. Zeldovich, ZhETF Pis'ma Red. 4 (1966) 174; 
R. Cowsik and J. McClelland, Phys. Rev. Lett. 29 (1972) 669. 

[3] A. D. Dolgov, S. H. Hansen, G. Raffelt and D. V. Semikoz, Nucl. Phys. B 580 
(2000) 331 ||hep-ph/0002223H ; 

A. D. Dolgov, S. H. Hansen, G. Raffelt and D. V. Semikoz, Nucl. Phys. B 590 
(2000) 562 ||hep-ph/0008T38H . 



[4] NOMAD Collaboration, |frep-ex/0101041 1 



[5] K.A. Olive and M.S. Turner, Phys. Rev. D25 (1982) 213. 

[6] S. Dodelson and L. M. Widrow, Phys. Rev. Lett. 72 (1994) 17. 

[7] S. Colombi, S. Dodelson, and L.M. Widrow, Astrophys. J. 458 (1996) 1. 



11 



[9 
[10 



[11 
[12 

[13 

[14 

[15 

[16 

[i7; 

[18 
[19 
[20 

[21 



R. A. Malaney, G. D. Starkman and L. Widrow, Phys. Rev. D52 (1995) 5480. 

G. Dvali and Y. Nir, JHEP 9810 (1998) 014. 

Z.G. Berezhiani and R.N. Mohapatra Phys. Rev. D52 (1995) 6607; 

Z.G. Berezhiani, A.D. Dolgov, and R.N. Mohapatra, Phys.Lett. B375 (1996) 26; 

Z.G. Berezhiani, Acta Phys.Polon. B27 (1996) 1503. 



J. Sommer-Larsen and A.D. Dolgov, ||astro-ph/9912166 
S. Hannestad, |astro-ph/ 000845 1 



X. Shi and G. M. Fuller, Phys. Rev. Lett. 82 (1999) 2832. 
M. Drees and D. Wright, ||hep-ph/ 0006274] . 



G. Fuller, "Comment on "New constraints on cool Dark Matter"", to appear. 

R. Barbieri and A.D. Dolgov, Phys. Lett. B237 (1990) 440; 
K. Kainulainen, Phys. Lett. B244 (1990) 191. 

K. Enqvist, K. Kainulainen, and M.Thomson, Nucl. Phys. B373 (1992) 498; 
A.D. Dolgov, ||hep-ph/0006T03H . 



D. Ndtzold and G. Raffelt, Nucl. Phys. B307 (1988) 924. 

V. Barger, R. J. Phillips and S. Sarkar, Phys. Lett. B352 (1995) 365. 

F.W. Stecker, Phys. Rev. Lett. 45 (1980) 1460. 

R. Kimble, S. Boyer, and P. Jacobsen, Phys. Rev. Lett. 46 (1981) 80. 



[22] E. W. Kolb and M. S. Turner, "The Early Universe," Redwood City, USA: Addison- 
Wesley (1990). 

[23] M. T. Ressell and M. S. Turner, FERMILAB-PUB-89/214-A. 

12 



[24] K. Kainulainen, J. Maalampi, and J.T. Peltoniemi, Nucl. Phys. 358 (1991) 435; 
G. Raffelt and G. Sigl, Astropart. Phys. 1 (1993) 165. 

[25] G. G. Raffelt, Chicago, USA: Univ. Pr. (1996) 664 P- 

Lacking a better synonym than lukewarm dark matter, we decide to use the standard 



name sterile neutrinos (see [11, 13 1 for great names like warmons or coolDM). 



[27] P. Colin, V. Avila-Reese and O. Valenzuela, Astrophys. J. 542 (2000) 622. 

[28] A. Kusenko and G. Segre, Phys. Lett. B396 (1997) 197; 
A. Kusenko, [|astro-ph/9903167| . 



13 



